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^ Abstract 

^ We report comprehensive high-accuracy molecular dynamics simulations using the ReaxFF 

^ forcefield to explore the structural changes that occur as Au nanowires are elongated, estab- 

lishing trends as a function of both temperature and nanowire diameter. Our simulations and 
IL* subsequent quantitative structural analysis reveal that polytetrahedral structures (e.g. icosa- 

^ hedra) form within the "amorphous" neck regions, most prominently for systems with small 

^ diameter at high temperature. We demonstrate that the formation of polytetrahedra diminishes 

the conductance quantization as compared to systems without this structural motif. We demon- 
strate that use of the ReaxFF forcefield, fitted to high-accuracy first principles calculations of 
Au, combines the accuracy of quantum calculations with the speed of semi-empirical methods. 
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The observation of conductance quantization in atomic-scale junctionsl^^ has fueled much in- 
terest in understanding their formation in mechanically deformed metallic nanowires (NWs),'^^ 
in large part due to potential applications in nano- and molecular-electronic devices. To this extent, 
quantifying the structural changes that occur as a NW elongates is important for understanding how 
and why parameters such as initial crystal structure,!^ temperature,!^ and rate of elongation'^^'^ 
affect the material properties (e.g., conductance,!^^^ mechanical stability,'^^'^!] qj. ^j^^ opening of 
a bandgapl^. Many studies have been performed that focus on identifying changes to the crys- 
talline structure, such as slipping and reorientation of the crystalline lattice,^ ^*^"^ * ^^*^"^ * and under 
what conditions single-atom wide monatomic chains are likely to form.f^EIIinini However, to date, 
there has only been limited emphasis on quantifying the formation and structure of the "amor- 
phous" or, more generally, non-crystalline domains that often appear within NWs.*^^ ^*^^ * ^^ l The 
formation of these non-crystalline structures is likely to result in significant changes to the proper- 
ties of the NWs. For example, amorphous structures separating otherwise crystalline NW domains 
may exhibit unique dynamical properties akin to glass-forming liquids, similar to atoms at grain 
boundaries in bulk metals. Moreover, non-crystalline structures may still demonstrate distinct 
local ordering, such as the distorted icosahedral structures predicted to form in the necks of Na 
NWs^^ and the recent synthesis of bi-metallic NWs with local icosahedral ordering. This is of 
particular consequence as thiolated icosahedral Au NWs are predicted to behave as semiconductors 
under certain charge states^ and non-crystalline NW neck structures have recently been demon- 
strated to produce gradual changes in the conductance, as opposed the more typical quantization, 
i.e., "steps" and "plateaus".*^ 

In this paper, we focus on identifying the amorphous/non-crystalline structures that form within 
simulated Au NWs undergoing tensile elongation. Simulation is well suited to study this problem 
as it provides the full spatial and temporal coordinates of the atoms in the NW, allowing the struc- 
ture to be quantified by the use of order parameters, such as those recently adapted from the field 
of shape matching .'^SHSl Much of the previous simulation work dealing with the elongation of Au 
NWs has relied on either quantum or, more commonly, semi-empirical 
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methods,!^ ^ ! ^ ! ^^*^^ * ^^*^^*^^ * ^^*^^ * requiring a compromise between the accuracy of the interatomic in- 
teractions and the accessible system size and timescale. Quantum mechanical methods, such as 
density functional theory (DFT), are widely considered to provide a highly accurate description of 
the interactions between Au atoms, in particular, capturing the stability of 2-d planar geometries 
for low coordination clusters .l^^l^ However, the high computational cost of quantum mechanical 
methods often limits both the total number of independent statepoints that can be efficiently consid- 
ered and the total size of the NW (typically ^100 atoms or less'^^'^, both of which may inadver- 
tently bias the results. Semi-empirical methods, such as the second-moment approximation of the 
tight-binding (TB-SMA) scheme,'^^! are many orders of magnitude faster than quantum mechanical 
methods, allowing for the efficient exploration of large numbers of statepoints and systems sizes 
more representative of experiment J^ * ^^*^^*^^ While methods such as TB-SMA have been shown to 
qualitatively match the behavior of experiment in side-by- side studies,*^ semi-emperical methods 
are often fitted to bulk systems, and thus may produce results that are quantitatively inconsistent 
with theory or experiment,'^ especially for low coordination atomic-scale contacts .'^^MQI partic- 
ular, the minimal energy 2-d planar structures of Au are not typically stable for semi-empirical 
methods which may incorrectly bias NW simulations towards forming 3-d isomers. 

In order to bridge the gap between semi-empirical and quantum mechanical methods, we use 
the ReaxFF reactive forcefield*^ which is designed to provide accuracy approaching quantum me- 
chanical methods, but at a substantially reduced computational cost. ReaxFF forcefields are de- 
rived by fitting analytical functions to DFT bonding curves. '^^'^ The parameters derived for Au 
by Keith, et al.^ have been demonstrated to satisfactorily capture the expected behavior of both 
bulk and low coordination states.f^In particular, this parameterization has been shown to closely 
reproduce the energy and stability of low and minimal energy clusters of Au.l^ In this work, we 
show that ReaxFF is a substantial improvement over TB-SMA for the simulation of NW structures 
by comparing the absolute and relative energy scaling to DFT calculations. We perform com- 
prehensive molecular dynamics (MD) simulations'^ with ReaxFF to investigate the behavior and 
properties of [100] orientated FCC Au NWs undergoing tensile elongation in vacuum, which we 



3 



also compare to simulations performed using TB-SMA. We report, for the first time, the formation 
of well-ordered polytetrahedral local structures (e.g., full and partially coordinated icosahedral^ 
within the non-crystalline domains of the NWs. We quantify the formation of these structures by 
using the Ryi^ shape matching 

method™Ii5] 

based on spherical harmonics,!^ which we use to 
construct trends as a function of temperature and NW diameter. We further calculate the zero- 
bias conductance of our simulation trajectories to demonstrate the impact that the formation of 
polytetrahedra has on the conductance behavior. 



Results 

Validation of ReaxFF for Au NWs 

We first assess the energetic description of Au NWs predicted by ReaxFF by comparing to DFT and 
TB-SMA. We report the energy for [100] oriented periodic FCC NWs with diameters, D=l.l, 1.5, 
and 1.9 nm, as well as an elongation sequence of a periodic 0.85 nm NW. The elongation sequence 
is representative of the low coordination structures typically found in the necks of mechanically 
deformed NWs. We perform DFT calculations that make use the generalized gradient approxima- 
tion by Perdew, Burke, Ernzerhof (PBE-GGA) and additionally include spin-orbital coupling (see 
Methods). All reported DFT energy values are calculated using structures that have been fully and 
self-consistently energy-minimized. 



In [figure] [!][]! a. we plot the energy calculated with DFT vs. the energy calculated with 



ReaxFF. In this plotting scheme, an ideal match between DFT and ReaxFF occurs if the data point 



falls along the line y = x. In [[figure] [!][]! ^ , filled symbols represent the energy calculated us- 
ing the coordinates energy-minimized with DFT (i.e., when calculating the ReaxFF energy, atoms 
are static and arranged according to the configuration generated with DFT). The DFT energy- 



minimized configurations for the elongation sequence are shown in Figures [[figure] [ 1 ] []T] :-g. We 



observe that ReaxFF is capable of accurately reproducing the per particle energy predicted by 
DFT, where the magnitude of the difference is less than 0.05 eV per particle and all data points 
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Figure 1: (a) DFT energy vs. ReaxFF energy, (b) DFT energy vs. TB-SMA energy. In both (a-b) 
filled symbols correspond to the energy calculated of structures that have been energy-minimized 
with DFT, open symbols correspond to energy of structures energy-minimized with each of the 
methods (e.g., DFT energy-minimized and ReaxFF energy-minimized). Snapshots of the energy- 
minimized configurations of the elongation sequence for: (c-g) DFT, (h-1) ReaxFF, and (m-q) 
TB-SMA. 



fall along the line y =x. In [[figure] [!][]! ^ , data points with open symbols correspond to energy- 
minimized structures predicted by the two different treatments (i.e., the DFT energy corresponds 
to the structure energy-minimized with DFT and the ReaxFF energy corresponds to the structure 
energy-minimized with ReaxFF). The ReaxFF structures are generated starting from the DFT op- 
timized structures. The energy-minimized elongation sequence for ReaxFF is shown in Figures 
[figure] [!][]! | i-l, where we observe minor changes as compared to the DFT structures in Figures 



[figure] [!][]! | :-g. Note, the preferred interatomic spacings will vary slightly between methods, thus 
we use a two stage energy minimization for ReaxFF, first allowing the box volume to change to 
minimize pressure (which scales atom positions), then allowing atom positions to change in order 
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to minimize the potential energy. For minimized configurations, ReaxFF predicts a lower energy 
structure than DFT for all statepoints. However the values are still close, where the magnitude of 
the deviation varies from 0.09 to 0.20 eV/particle. 



[figure] [!][]! 3 is constructed in the same manner as [figure] [!][]! a, but using TB-SMA rather 



than ReaxFF. Again, filled symbols represent the energy calculated using the DFT energy-minimized 
structures. TB-SMA fails to closely reproduce the energy, predicting lower values than DFT for all 
structures. Here the magnitude of the differences ranges from 0.17 to 0.41 eV per particle. Open 
symbols again represent the energy of structures minimized with the two methods. We use a sim- 
ilar procedure for TB-SMA as we did for ReaxFF. The energy minimized structures of TB-SMA 
are shown in | [figure] [!][]! ) n-q, where we observe only minor differences when compared to either 
DFT or ReaxFF structures. The most noticeable difference is observed in [[figure] [!][]! where 
the left hand side of the NW compresses into a more locally dense state. For energy-minimized 
structures, the energy predicted by TB-SMA is shifted to even more negative values, where now 
the magnitude of the energy difference with DFT ranges from 0.21 to 0.47 eV per particle. 

In addition to comparing the absolute energy, we examine the relative change in energy ranging 
from the highest coordination state (D=1.9nm NW) to the lowest coordination state (i.e., the final 



elongated structure, shown in Figures [[figure] [ 1 ] []Tg , 1 and q). This can be assessed by examining 
the slope of the data points shown in Figures [[figure] [ 1 ] [JT^ and b; a slope of unity provides an ideal 



energy change with respect to DFT. For the ReaxFF calculations in [[figure] [!][]! ^ , we observe a 
slope 0.99 for the static structure (filled symbols) and 1.11 for the energy-minimized structures 
(open symbols), using linear regression. Thus, even though the absolute energies are slightly 
different for the energy-minimized states (open symbols), the relative change predicted by ReaxFF 
is in very good agreement with DFT. For the TB-SMA calculations in |[figure][l][]l) 3, we observe 
a slope of 1.63 for static calculations (filled symbols) and a slope of 1.81 for energy-minimized 
structures (open symbols). The strong deviation from unity indicates that TB-SMA fails to properly 
predict the energy loss associated with a reduction in average coordination number, which may 
strongly influence the structural evolution of the NW. In previous work comparing common semi- 
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empirical methods for the elongation of finite Au NWs, TB-SMA was demonstrated to provide the 
best agreement with DFT in terms of the relative energy change.^i^l Thus it appears that ReaxFF is 
an improvement over TB-SMA, providing a better description of both the relative energy change 
and absolute energy for Au NWs. 

The formation of polytetrahedral structures 
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Figure 2: Representative ReaxFF simulation snapshots of NWs at T=298K elongated at 0.1 m/s 
for (a) D=l.lnm NW elongated by 12.2 A and (b) D=1.5nm NW elongated by 21.9 A . In a 
and b, atoms are rendered at ^ 75% of their true size with individual polytetrahedral structures 
highlighted within the NW and shown at the right, (c) Rendering of a Z12 structure, i.e., an 
icosahedron, with an individual tetrahedral subunit highlighted. Partial Z12 clusters are shown as a 
function of coordination number, cn, varying from fully coordinated, cn=12, to cn=6. Note, cn=6 
shows two different partial clusters commonly observed. 



We perform simulations of cylindrical Au NWs of length 3.6 nm, with diameters, D=l.l, 1.5 
and 1.9 nm (177, 317, and 589 atoms, respectively) for a variety of statepoints, calculating 10 
independent simulations for each parameter combination (see Methods for more details). Here 
we focus on ReaxFF-based simulations but also include TB-SMA-based simulations for compar- 
ison. In these ReaxFF simulations we observe that, as a typical Au NW elongates, a small region 



7 



begins to "neck," becoming narrower than the original diameter. The neck forms in response 
to the apphed tensile load, localizing many of the structural changes to this neck region, allow- 
ing the bulk of the NW to relax its configuration. Within the necks of our simulated NWs, we 
often observe the formation of non-crystalline structures with distinct ring- and crescent-like ge- 
ometries. I [figure] [2] []2^ shows a representative ReaxFF simulation snapshot of a D=l.l nm NW 
at T=298K, shortly after neck formation; an individual ring-like structure is highlighted in the 



snapshot with an enlarged version, extracted from the system, shown to the right, [[figure] [2] []2\ ) 
shows a representative snapshot of a D=1.5nm NW where two crescent-like structures are high- 
lighted and extracted. These ring- and crescent-like structures are polytetrahedral in nature, where 
the local structure of the atoms can be decomposed solely into individual tetrahedral subunits.!^ 
Specifically, the highlighted structure in [[figure] [2] []2^ is best classified as an icosahedron. An 
icosahedron is constructed of 13 particles arranged as two splayed pentagonal pyramids (com- 
posed of 20 total tetrahedral subunits), demonstrating 5-fold symmetry;S51an ideal icosahedron is 
shown in [[figure] [2] []2[ :. The two highlighted clusters in [[figure] [2][]2} 3 are partially coordinated 
icosahedral clusters, that maintain the same basic geometry of the full cluster, but with X number 
of particles removed (i.e., several tetrahedral subunits removed) .f^Sl Ideal partial clusters of various 



coordination are shown in [[figure] [2] []2| :. TB-SMA simulations also demonstrate polytetrahedral 
structures within the necks of the elongated NWs, similar to Figures [[figure] [2] []2^ and b. 

To quantify the structures formed in our simulations and establish trends, we use the Ryi^n 
method to determine which atoms are in polytetrahedral local environments (see Methods). We 
use this quantitative information to construct three metrics: (1) the maximum number polyte- 
trahedra that occur concurrently during the course of a simulation trajectory, (2) the number of 
polytetrahedra per timestep in a simulation (calculated by normalizing the total number of polyte- 
trahedra observed by the total number of statepoints considered), and (3) the coordination number 
(cn) histogram of polytetrahedral atoms. 

We first report the behavior of thin NWs with D=l . 1 nm elongated at a rate of 1 m/s for T={ 10, 
100, 200, 298, 400, 500} K. Figures pgure] [3] []3^ and b plot the maximum number of polyte- 
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Figure 3: Quantification of polytetrahedral structures, (a) The maximum number of polytetrahedral 
clusters observed concurrently in a trajectory and (b) the number of polytetrahedra per timestep, 
both for D=l . Inm and shown as a function of temperature for NWs pulled at 1 m/s, calculated from 
120 total independent simulations, (c) The maximum number of polytetrahedral clusters observed 
concurrently in a trajectory and (d) the number of polytetrahedra per timestep, both at T=298K 
and shown as a function of diameter, calculated from 400 total independent simulations. Data for 
(c-d) is presented as the aggregate of simulations conducted from 0.1 to 5 m/s, as a strong rate 
dependence was not observed. All data points correspond to the median and error bars correspond 
to 1st and 3rd quartiles. The legends shown at the top are applicable to all plots. 



trahedra and the polytetrahedra per timestep, respectively. In all cases, we report median values, 
where error bars correspond to the first and third quartiles, since the data is not necessarily well 
described by a Gaussian distribution. ReaxFF simulations demonstrate a steady increase in both 
metrics as temperature is increased, i.e., both the maximum size of the polytetrahedral domains 
and the likelihood of finding polytetrahedra increases as T increases. The polytetrahedral motif is 
not prominent for ReaxFF systems at low temperature, i.e., T < 200K. While the behavior of the 
ReaxFF and TB-SMA potentials converge for high temperature, TB-SMA demonstrates markedly 
different trends and behavior. Overall, we observe that TB-SMA predicts a significantly larger 
value of both the maximum number of polytetrahedra and polytetrahedra per timestep than the 
ReaxFF simulations. TB-SMA predicts two regimes with an approximate crossover of T=200K. 
TB-SMA systems at T=10 and lOOK tend to maintain a predominantly FCC structure, whereas 
higher temperature states transition into A3 core-shell helical structures'^ upon moderate amounts 
of elongation. For example, at T=298K, 70% of the TB-SMA simulations transform into predom- 
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inantly A3 structured NWs by the time the NW has been elongated by 5A. The A3 structure is 
locally composed of tetrahedral subunits that closely match the polytetrahedral structures in our 
Ryiyyi refcrencc library. However, we note that even for T<200K, where A3 structures are not found 
with either potential, ReaxFF still forms fewer polytetrahedra than TB-SMA, by a factor of ^10 for 
both metrics. The predominance of crystalline structures at low temperature predicted by ReaxFF 
is more consistent with results of zero temperature DFT calculations performed for the elongation 
of [110] and [1 1 1] Au NWs.™] 

Figures [[figure] [3] []3| : and d plot the median value of the maximum number of polytetrahe- 
dra and polytetrahedra per timestep, respectively, for NWs as a function of NW diameter; these 
simulations are performed at T=298K, as this statepoint demonstrates appreciable numbers of poly- 
tetrahedra for both potentials and temperatures in this vicinity are common in experiment.'^ The 



values reported in Figures, [[figure] [3] []3| : and d correspond to the aggregate of simulations per- 
formed for rates 0.1 to 5 m/s, as we did not observe strong variability with rate of elongation. 
We observe that the median value of both metrics decreases with increasing diameter for both 
ReaxFF and TB-SMA. That is, the likelihood of forming polytetrahedra decreases with increas- 
ing NW diameter. In all cases, ReaxFF-based simulations predict smaller median values for the 
two metrics than TB-SMA, with the largest differences occurring for the D=l.lnm NWs. The 
differences we observe between ReaxFF and TB-SMA are likely related to the energy scaling we 
previously explored, where TB-SMA overpredicts both the absolute magnitude and relative energy 
of lower-coordination states. 



[figure] [4][]4 plots histograms of the coordination number of atoms identified as having poly- 



tetrahedral local ordering for NWs at T=298K (corresponding to the data presented in Figures 



[figure] [3] [] 3: and d). Histograms of ReaxFF simulations demonstrate a diameter dependence; 



D=l.lnm NWs show two key peaks at cn=6 and 8, but the relative fraction of cn=6 strongly in- 
creases with increasing diameter. For D=1.9 nm NWs simulated with ReaxFF, ^ 80 % of the 
observed polytetrahedral clusters have cn=6. Histograms generated from TB-SMA simulations 
demonstrate two clear peaks at cn=6 and 8, encompassing ^ 80 % of the total polytetrahedral 
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Figure 4: Histograms of the coordination number (cn) of polytetrahedral structures for NWs with 
(a) D=l.lnm, (b) D=1.5nm, and (c) D=1.9nm. All data is at T=298K and is presented as the 
aggregate of simulations conducted from 0.1 to 5 m/s. The legend at top is applicable to all plots. 



structures, with much less dependence on diameter. 



Zero-bias conductance 

We calculate the zero-bias conductance of several representative trajectories generated with ReaxFF 
and TB-SMA in order to assess the impact of polytetrahedra, i.e., is there a characteristic conduc- 
tance "fingerprint" associated with their formation. As described in the Methods section, we use a 
combined DFT-non-equilibrium Green's function technique to calculate the conductance,^^ which 
we report in terms of Go = 2e^/h, where e is the charge of an electron and h is Planck's con- 
stant. In experiment, a typical conductance trace demonstrates a "staircase" behavior, with abrupt 
plateaus and drops as a function of elongation. The abrupt changes in the conductance coincide 
with structural rearrangements of the atoms, such as thinning of the NW neck. For very small 
cross-sections, the steps often occur at roughly integer multiples (m) of Gq. It is important to note 
that histograms of conductance tend to demonstrate a substantial spread around each of these inte- 



11 



ger peaks'^ and, as a result, non-integer values are also regularly observed in experiment. Recent 
work has correlated the formation of non-integer plateaus to the difference in dimensionality of the 
local structures within the NW neck (i.e., the difference between 3-, 2-, and 1-dimensional atomic 
arrangements) ^ 
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Figure 5: (a) Conductance and (b) number of polytetrahedra as a function of elongation of a 
ReaxFF NW at T=298K. (c) Conductance and (d) number of polytetrahedra of a ReaxFF NW at 
T=10K. (e) Conductance and (f) number of polytetrahedra of a TB-SMA NW at T=10K. All data 
corresponds to NWs with D=l.lnm elongated at 1 m/s. Plateaus are highlighted with horizontal 
dotted lines. Simulation snapshots for configurations with G^SGq are inset in each figure, where 
c and e each highlight a polytetrahedral structure in the neck (red). 



[figure] [5][]5| plots the conductance and number of polytetrahedra as a function of elongation 
for trajectories of D=l.lnm NWs elongated at 1 m/s, generated with both ReaxFF and TB-SMA. 
We first consider a ReaxFF trajectory at T=298K that does not demonstrate any significant polyte- 
trahedral ordering, as shown in Figures [[figure] [5][]5^ -b. We observe distinct step-wise behavior 
with conductance plateaus at roughly integer multiples of Go (m ^ 1,2, 3 and 4 Gq, highlighted in 



[figure] [5][]5^ with horizontal lines). Additionally, the length of these plateaus range from 1.5-3.5 

O OA 

A which matches the approximate plateau length of 1-4A seen in experiments at T=4K and 



T=298K.^ A snapshot of the NW at ^3 Go is inset in [figure] [5][]5a. showing a crystalline neck 
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structure. 



Figures [[figure] [5] []5p -d plot data for a NW trajectory generated with ReaxFF at T=10K. As 
before, we observe step wise behavior in the conductance. In contrast to the previous trajectory, 
we now observe the formation of an appreciable number of polytetrahedral structures, as shown 
[figure] [5][]5| i. The formation of polytetrahedra appears to have a distinct impact on the con- 



in 



ductance. For example, while we observe plateaus at m ^3 and 1 Go in [[figure] [5][]5| :, there is a 
gradual change in conductance between these two steps rather than a plateau at 2 Go . The onset 
of the formation of polytetrahedral structures at ^5A correlates with the start of the conductance 
plateau at 3Go. Visually, we observe that at the start of the gradual change (i.e, 3Go), the narrowest 
region of the neck is spanned by a polytetrahedral structure ^ 3 atoms wide; a shapshot is inset in 



[figure] [5] []5 



We further assess the impact of polytetrahedra by examining a trajectory at T=10K generated 
with TB-SMA, as shown in Figures pgure] [5] []5^ -f . Recall that TB-SMA over-predicts the ten- 
dency to form polytetrahedra as compared to ReaxFF, making it ideal for assessing their impact. 
We still observe step-like behavior, but now only see plateaus at m ~ 4 and 1 Go. Again, we ob- 
serve a rather gradual transition from 3Go to IGo. While the onset of polytetrahedra appears at 2.5 



A of elongation, the most significant peak occurs at ^5 A (see [figure] [5] [] 5 which correlates 



to the start of the gradual change in conductance. As was noted in |[figure] [5][]5| :, we also observe 
that at ^3Go (i.e., the start of the gradual change), the narrowest region of the neck is spanned 
by a polytetrahedral structure ^ 3 atom diameters across (a shapshot is inset in [[figure] [5][]5| :). 
Thus it appears that the gradual conductance changes we observe are associated with the forma- 
tion of polytetrahedra that span the cross-section of the neck. This is likely related to different 
modes of mechanical deformation exhibited by tetrahedral structures verses crystals; the mechan- 
ical behavior of related core-shell NWs as compared to crystals has previously been explored in 
RefP 

Our results regarding the impact of polytetrahedral neck structures on the conductance is sup- 
ported by recent simulations of Tavazza, et al that used DFT to investigate the elongation of [1 10] 
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and [1 1 1] Au NWs at zero temperature.'^ In that work, nearly identical conductance behavior was 
observed for systems that demonstrated non-crystalline neck structures, i.e., gradual changes from 
3Go to IGq.^^ Similar conductance traces were also observed in experiments at T=4K.^^ Addition- 
ally, it was noted that systems exhibiting this gradual conductance change occurred in the minority 
of simulations and experiments, suggesting that non-crystalline structures are rare at these very 



low temperatures. '22 As previously shown in Figures [figure] [2] []2 i and b, our ReaxFF simula- 
tions also predict that the non-crystalline polytetrahedral structures are a rare structural motif at 
low temperature. 



Discussion 

The formation of polytetrahedral structures within the necks of the elongated Au NWs may be 
surprising as icosahedra are not predicted to be the lowest energy structures. I^^^ Generally speak- 
ing, planar structures are predicted to be the minimal energy states for small, isolated clusters 
of Au.'^^'^Ql However, this does not imply that icosahedra, or other low, but not minimal energy 
structures are inherently unstable. To assess the stability of the polytetrahedral structures, we 
perform conjugate gradient energy minimization using DFT of the isolated icosahedral clusters 



shown in [[figure] [2] []2| :. Starting from both ideal structures and those generated with ReaxFF, 
we find that clusters with cn=6, 8, 9, and 12 maintain their original geometry at the end of the 
minimization, suggesting that these structures are locally stable. We also note that systems that 
start as partial icosahedra with cn=10 and 11 relax their structures into more "flattened" states 
that, in 3 out of 4 calculations, contain a distinct polytetrahedral pentagonal di-pyramind (i.e., a 
cn=6 polytetrahedral structure). Recall that the cn=6 structure is the most common polytetrahe- 
dral configuration found in our simulations. Furthermore, polytetrahedral structures of Au are 
well reported in experiment and theory, including monolayer protected icosahedral clusters with 
cn=l 2,1221 phosphone-thiolate-protected partial-polytetrahedral Au structures,!^ and bi-icosahedral 
clusters of both face-sharinglSSI and vertex-sharing icosahedra.l^Sl 
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It is important to note that the polytetrahedral structures we observe in our simulations are not 
isolated clusters, but rather occur within the context of a larger, mechanically deformed NW. The 
average coordination number of atoms in an isolated cluster will be lower than an equivalent config- 
uration that occurs within a NW. Additionally, we must also consider that these structures form as a 
result of mechanically deforming a NW via an applied tensile load. Previous experiments of NWs 
pulled at a constant rate demonstrate a saw-tooth like behavior in the mechanical force on the NW 
as a function of elongation and thus the local forces of the atoms may be far from equilibrium 
and constantly changing. The original structure of the NW will also dictate how it can respond to 
tensile loads.'^'^For example, it has been demonstrated that the likelihood of forming monatomic 
Au chains in elongated Au NWs is influenced by the original orientation of the crystalline lattice 
(i.e., [Ill] vs. [110] vs. [100]). 13 Focusing only on minimal potential energy structural motifs also 
ignores the importance of thermal effects. Recall that in Figures [[figure] [3] []3^ -b we observed that 
the formation of polytetrahedral structures is very unlikely at low temperature, but increases with 
temperature. Previous work from our group proposed a universal energy release mechanism that 
outlines the strong role that temperature can play in determining how a NW responds to the applied 
tensile load.^ 

Additionally, Au is ductile, which manifests itself as a tendency for the NW to remain in a 
cylindrical, connected structure by thinning within the neck region. As the neck of the NW thins, 
polytetrahedral structures may become the optimal configurations, since they are able to span nar- 
row, cylindrical-like geometries by forming face-sharing polytetrahedral or vertex- 
sharing structures, was demonstrated for particles with van der Waals interactions, where 
a transition was observed from FCC/HCP structures to polytetrahedral structures upon decreasing 
the diameter of the confining cylinder. Similar behavior was also demonstrated for Au, where 
tetrahedral helices were predicted to be the optimal structure for small diameter NWs whereas 
FCC structures were optimal for large diameter NWs.*^ Recent experimental work has demon- 
strated the synthesis of thin Au/Ag NWs that are locally constructed of face- sharing icosahedra,'^^ 
additionally demonstrating this packing motif for metallic NWs. Furthermore, the formation of 
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polytetrahedral structures within the necks of elongating NWs is not without precedent. Previous 
DFT calculations of the mechanical deformation of [100] oriented BCC Na NWs demonstrated the 
formation of distorted full and partial icosahedra within the NW neck. This may suggest that the 
formation of polytetrahedra is a more universal relaxation mechanism for mechanically deformed 
NWs. 

The formation of polytetrahedral structures may have important implications regarding the 



electronic properties of the NW. As we saw in [[figure] [5][]5[ the formation of polytetrahedra di- 
minishes the quantization of the conductance traces. Combined with the increased likelihood of 
forming polytetrahedra as temperature increases, this may be an important to consider when en- 
gineering nanoelectronic devices designed to operate at, e.g., room temperature. The impact of 
polytetrahedra may be even more substantial. The work of Jiang, et al predicts that thiolated 
vertex- sharing icosahedral Au NWs behave as a semiconductor for the -2 and +2 charge states.l^Sl 
This suggest that Au NWs elongated in the presence of a thiol solvent might, under certain condi- 
tions, change from behaving as metals to semiconductors as a result of the formation of polytetra- 
hedral neck structures. 



Conclusion 

In this work, we showed good agreement between ReaxFF and DFT in terms of the absolute energy 
and relative energy scaling. We also showed that ReaxFF provides a better energetic description 
than the semi-empirical TB-SMA potential, which has commonly been used to study Au NW 
elongation. We further demonstrated the formation of polytetrahedral local structures within the 
necks of elongated Au NWs simulated with ReaxFF. We observed that the formation polytetrahedra 
is not a strong motif of structural relaxation for low temperature systems, but the likelihood of 
forming polytetrahedra increases with temperature. We also found that the likelihood of forming 
polytetrahedra increases as NW diameter decreases. Polytetrahedral order was also observed for 
simulations performed with the TB-SMA potential, however, TB-SMA appears to significantly 
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over-predict their formation as compared to ReaxFF, especially for small diameter NWs and low 
temperature. We find that ReaxFF-based simulations produce trajectories that demonstrate clear 
step-wise behavior in conductance that closely matches experiment. We also observe that the 
formation of polytetrahedra (i.e. non-crystalline structures) within the NW necks dampens the 
quantization behavior of the conductance, closely matching the observations in Ref The results 
presented here suggest that the formation of polytetrahedral structures within the necks of the 
[100] orientated Au NWs is a ubiquitous and important mechanism for relieving the internal stress 
brought about by elongation, in particular, for small diameter NWs at high temperature. 
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Methods 

DFT Calculations: We perform DFT calculations using the generalized gradient approximation of 
Perdew, Burke, Ernzerhof (PBE-GGA)!^ within the plane wave projector augmented wave formal- 
ism (PAW).'^^'^ Our calculations were performed at the gamma point with a 230eV plane- wave 
energy cutoff. The PAW potentials included the outermost d and s electrons in the valence, giving 
1 1 electrons per atom. Significantly, we fully included spin-orbit coupling in all our calculations 
of the energy and forces. All reported DFT energy values are calculated using structures that have 
been fully and self-consistently energy-minimized. 

Simulation Methodology: Our basic simulation procedure closely mimics the "stretch and re- 
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lax" procedure of Pu, et al. ^ we rigidly fix the last two rows of each end of the NW creating 
"grip" atoms at each end; we allow the core of the system to relax for 100 ps; after relaxation, 
we elongate the NW by displacing the grip atoms at one end by increments of 0.05 A , allowing 
the system to relax for a prescribed time between displacements, repeating until the NW breaks. 
We explore elongation rates from 5 m/s (1 ps relaxation between displacements) to 0.1 m/s (50 
ps). This range encompasses typical elongation rates used to draw metallic whiskers with STM or 
Simulations are conducted between T=10 K and 500 K; common experimental conditions 
are T=4-10 K'^'^and SOOK.^^For each statepoint we perform 10 independent simulations to 
obtain sufficient statistics. 

ReaxFF simulations use the fitting parameters derived by Keith, et al.l^ This particular fitting 
was chosen because it was parameterized using PBE-GGA DFT, which is well suited for systems 
with surfaces (i.e., NWs), and the elastic moduli of the fitting closely match experiment.'^ TB- 
SMA simulations are performed using fitting parameters derived by Cleri and Rosato.^^ TB-SMA 
is used as previous comparisons with DFT showed it provides a better match in terms of relative en- 
ergy change than other common semi-empirical potentials.'^ Both ReaxFF and TB-SMA simula- 
tions are performed using the LAMMPS simulation package!^ (extended to incorporate TB-SMA) 
with the Nose-Hoover thermostat for systems without periodic boundary effects. As is common, 
TB-SMA simulations use a timestep of 2 femtoseconds,'^'^^'^^'^ whereas a smaller timestep of 
0.5 femtoseconds is needed for ReaxFF. To simulate the same timescale, ReaxFF requires more 
walltime than TB-SMA, by a factor of --50. 

Ryiyyi Method: To analyze the local configurations of atoms, we employ techniques borrowed 
from the field of shape matching }^^^^ We use the Ry^ scheme proposed in referencesj ^^'^^' which 
allows us to objectively identify structural trends in an automated fashion. This method relies 
on creating a rotationally invariant spherical harmonics'^ "fingerprint" of the first neighbor shell 
surrounding each atom. We compare the fingerprint of each local cluster to a library of reference 
configurations and determine the best match. Here, we quantify the match, M, using the Euclidian 
distance, normalized such that it spans from (worst match) to 1 (perfect match). By definition. 
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every cluster in this scheme has a "best match" even if that match is poor. To avoid misidenti- 
fication, we define a cutoff such that we ignore matches where M < 0.86. This provides results 
consistent with visual inspection; the choice of cutoff will impact the specific calculated values, 
however the overall trends are not strongly perturbed by this choice.^ Our reference library in- 
cludes full and partially coordinated clusters of type FCC, HCP, and icoshaedra, incorporating both 
ideal configurations and configurations relaxed using the ReaxFF and TB-SMA potentials, where 
we only consider clusters of 6 < cn < 12. 

Conductance Measurement: Conductance calculations are carried out in a combined DFT-non- 
equilibrium Green's function framework. A double-zeta numerical atomic orbital basis set is used 
(18 basis functions per Au atom) in the calculations. The zero-bias conductance is given by G = 
T{Ef) where T{E) is the transmission coefficient of the device and Ep is the Fermi energy. G is 
reported in units of Go = 2e^ /h, where e is the charge of an electron and h is Planck's constant. We 
validated our conductance measurements by investigating single point atom contacts, finding good 
agreement with previous theoretical calculations. '^^'^^J^ Further details of the implementation and 
its accuracy can be found in Ref.EH 
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Polytetrahedral structures form as Au nanowires undergo elongation. 
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